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We show that in attractive force fields a high speed particle with a small mass speeding through 
space, statistically loses energy by colliding softly and transversely with large masses that move 
slowly and randomly. Furthermore, we show that the opposite holds in repulsive force fields: the 
light particle statistically gains energy. We consider the large masses to be distant from one another 
such that a series of two-body scatterings would suffice for the statistical analysis. The difference 
of this process from equilibration is discussed. 

PACS numbers: 



Dynamical Description- In addition to the well known 
gravitational and Coulomb, nearly all other interactions 
in nature such as intermolecular forces and interaction of 
vortices in superconductors are central. In conservative 
fields, the central force on each particle can be derived 
from a potential function by F — — W(r) where r is 
the distance between the two bodies and V(r) = -p?; 
here a is the strength of the interaction depending on 
the parameters of the problem, k defines the range of the 
interaction [l|, 0| • Lastly, a < and a > correspond to 
attractive and repulsive force fields respectively [3j. 

In this paper we study a fast small mass passing 
through a dilute system of randomly moving central 
forces, where changes in the state of the small mass can 
be well approximated by a series of two-body scatterings. 
For the purposes of this section, it is best to think of the 
scatterings in the lab frame. Consider a scattering be- 
tween two interacting particles mi and m,2 where mi is 
much lighter (mi -C m-i) yet much faster than the other 
(v\ 3> V2) but nevertheless m\v\ <C miv\. For example 
one can visualize a small comet (mi) undergoing a small 
angle scattering in the gravitational field of the planet 
Jupiter (m.2). For the sake of concreteness take the po- 
tential to be attractive for now. In a typical scattering 
both masses are initially moving. The questions we are 
interested in investigating are: What statistically invari- 
ant features are shared by such scatterings in central force 
fields? Would many such small angle scatterings have a 
net effect on the energy of mi? 

While there is an active frontier of numerical work on 
many-body simulations [4] , there are still interesting sta- 
tistical inferences that can be derived from close analysis 
of two-body collisions. Further, such inferences in the 
cases where the collisions are head on have been exten- 
sively studied as in the Fermi acceleration mechanism 
In contrast the main contribution to the effect presented 
in this paper is from collisions that are transverse; i.e., 
the trajectory of the massive particle, roughly speaking, 
is perpendicular to that of the small particle. 

The organization is as following: 1) We give a heuris- 
tic description of the effect, 2) We provide a numerical 
verification 3) We prove the result by giving a mathe- 



matical derivation, 4) Lastly we discuss this effect in a 
larger statistical context and discuss its difference from 
the equilibration process. 

Let us consider two extreme cases that would convey 
the gist of what underlies this work. In the first case the 
massive particle, ni2, slowly and transversely veers away 
from the trajectory of mi that is speeding by. In the 
second case, m.2 slowly and transversely approaches the 
trajectory of mi (see Figure [3]-a). 

In the first case where m.2 is moving away, mi falls into 
the potential well of and so long as it is approaching 
the point of minimum distance it gains kinetic energy. 
After passing this point, mi starts climbing up the po- 
tential well and pays back the gained kinetic energy by 
restoring it into the potential energy of the two-body sys- 
tem. However, on the way out of the well it climbs a well 
that is effectively smaller than the one it fell into as m.2 
is on average farther away from it. Therefore, in the case 
that the big mass is moving transversely away from the 
particle, the small particle emerges with a gain in the 
kinetic energy i.e. ±m,i\v mi \ 2 _ oc < \mi\v mi \ 2 +00 . 

The exact opposite effect holds in the second case, 
where m.2 is moving towards mi. In this case, mi enters 
the potential well set up by m.2 and, as in the previous 
case, gains kinetic energy so long as it is approaching the 
minimum distance between the two masses. However, on 
the way out of the well it faces a more demanding climb 
as m-2 is on average closer to it and the potential well is 
steeper and deeper than before. Therefore, in the case 
that the big mass is transversely approaching the particle 
(Figure [lj, the small particle emerges with a loss in the 
kinetic energy i.e. ±mi\v mi ^ > 5TOi|i> m il+oo • 

The point however is that the two cases are not sym- 
metric. The decreasing of the magnitude of the force with 
distance breaks the symmetry between the two cases. This 
is shown in Figure [TJ In an attractive force field, mi 
has a greater loss (in magnitude) of energy when mi ap- 
proaches it than a gain (in magnitude) when m-2 recedes 
away from it. This asymmetry, deduced from dynamical 
principles, has consequences for the statistical mechanics 
of mi that we will explore. 

In repulsive force fields, a > 0, the potential in Fig- 
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Schematics of The Loss and The Gain in an An Attractive Potential 
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Figure 1: A plot of potentials of type V(r) — 7712 3> mi, 
originally sitting at r = 0, is the central force that mi scatters 
from. When 7712 moves to the right (approaches mi) by Sr, it 
will set up a potential that is AEi deeper for mi to climb out 
of than the one it fell into, whereby mi loses AEi of energy 
to the central force as the result of the scattering. Whereas, 
when 77i2 is moving to the left (away from mi) by the same 
amount, Sr, its potential is AE2 shallower and consequently 
mi gains AE2 of energy. Note that the nonlinear dependence 
of V(r), on r makes the magnitude of the energy transfer in 
the two cases asymmetric: |Ai?i| > | All- 
ure [T] flips about the horizontal axis. Therefore the phe- 
nomenology is the exact opposite. Namely, the particle 
has a larger gain than loss. 

So far we have described a purely dynamical phenom- 
ena where 777,2 collides softly and transversely with mi. 
Clearly if a very fast small particle (e.g., an electron) zips 
through a dilute soup of big masses (e.g, massive ions or 
stars in a galaxy) that randomly either approach it or 
move away from it in the manner described above, the 
small mass statistically loses (gains) energy to (from) the 
big masses when the force fields are attractive (repulsive). 

The fact that a high speed particle with low energy 
can lose or gain energy by scattering from larger masses 
with higher energies is not new. In the gravity assist or 
the slingshot effect, a satellite can be carefully steered so 
that, by elastically scattering from a planet or a star that 
is moving along or opposite to the direction of its asymp- 
totic motion (t — > 00), it acquires a boost in energy or 
brakes to save fuel The intuition in gravity assist, as 
in this work, is that to boost the energy of the satellite we 
have to throw it behind the approaching planet and vice 
versa Q- Slingshot effect, however, is not a net statistical 
effect resulting from an underlying random process. 

We proceed to examine a scenario where the large 
masses move in random directions. 

A Numerical Experiment- We simulate a process where 
Si = 5 x 1(T 8 , a = -0.13346, and k = 1. The initial 
conditions for mi and initial position of 7772 are taken 
to be the same in all the simulations (see Appendix U 
for the numerical details). To single out the effect of the 
direction of the motion of the large mass on the statistical 
energetics of mi, we take the velocity of 7772 to be: v^ 2 = 
V2 cos 9 and = v 2 sin 9, where t>2 is kept constant, 



Figure 2: The black line corresponds to the case that 7712 
is pinned down and mi, considered alone, undergoes a con- 
servative elastic scattering: v mi (+00) = v mi (— 00). The red 
line corresponds to the 6>— average of the asymptotic change 
in the kinetic energy of mi, with 7712 moving in twenty dif- 
ferent directions chosen uniformly on a circle. We find that 
AE = _E avg - E m2 Stationary = -6.441 x 10~ r . Error bars on 
the energies are of O(10~ 12 ) ; see the appendix. 

v\/v% ~ 10 and 9 is uniformly distributed on a circle 
9 = [0, jq, 2tt) . We find that mi on average loses 
energy to 7772 (Figure [3]). It is instructive to look at the 
points corresponding to the two extreme cases discussed 
in the previous section. When 7772 moves towards the 
trajectory of mi along the path of minimum approach 
(6 ~ 7T ), maximum loss is seen; furthermore, when 7712 
veers away from the path of minimum distance (9 ~ ^), 
the maximum gain is seen. Both are congruent with our 
expectations. In Figure [2j the ^-average of the change in 
the kinetic energy of mi is shown in red, which is seen to 
be less than the kinetic energy of mi when 7772 is pinned 
down (i.e. 7^2 = for all times). In the latter case, used 
as a calibration here, the asymptotic kinetic energy of mi 
is conserved (shown with a black line in Figure [5]). 

Analytical Derivation- We are considering an standard 
elastic collision between two particles [8|. Let Vi and 
V2 be the velocities of the two bodies in the lab frame; 
the velocities in the center of mass (C system denoted 
by the suffix 0) are related to the velocities in the lab 
by vio = m 2 v/(mi + m 2 ), v 2 o = —mxv/(m\ + m 2 ), 
where v = vi — \r 2 . In the C system the collision sim- 
ply rotates the velocities, which remain opposite in di- 
rection and unchanged in magnitude. If we denote by 
no the unit vector in the direction of the velocity of mi 
after the collision, then the velocities of the two par- 
ticles after the collision (distinguished by primes) are 
v' 10 = m2wn /(mi+m 2 ) and v 20 = -mivn /(mi +777 2 ). 
In order to obtain the final velocities in the lab frame, de- 
noted by v' x and v 2 , we add the velocity of the center of 
mass V = (mvi + 7772 V2) / (mi + 7712) to v^ and v 20 , 

v' x = m 2 wn /(mi + 777 2 ) + V, (1) 
v 2 = -777iTm /(mi + m 2 ) + V. (2) 
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Figure 3: Unbounded trajectory in a) an attractive force field 
and b) repulsive force field 



No further information about the collision can be ob- 
tained from the laws of conservation of momentum and 
energy. The direction of the unit vector no depends on 
the particular law of interaction and positions during the 
collision. We assume the massive particles are far enough 
from one another that a sequence of two-body scatter- 
ings would be adequate to describe the process [9]. The 
change in the kinetic energy of mi before and after a 



collision is AKE 



which reads 



a ( 

2 \ V 1 

AKE mi = (ivV. [nj - no] = fivV.n 



(3) 



where, a = mi . m2 is the reduced mass, n n and nt de- 
note the unit vectors pointing in the direction of motion 
of mi before and after the collision in the C frame and 
n = nj — n.Q is their difference. Before calculating Eq. [3] 
as a function of 9, once again, we look at the two special 
cases discussed above. First consider an attractive force 
field. Suppose V lies on the line of the minimum distance 
as shown in Figure [5^. Clearly if mi recedes away from 
m i then V points in the same direction as n and the 
dot product on the right hand side of [3] evaluates to be 
positive, whereas if V and n point in opposite directions 
the right hand side is negative. In the case of repulsion 
the signs would be the opposite (see Figure [3]d). 

What we now show is that the two cases are not sym- 
metric. That is the kinetic energy loss (gain) in the ap- 
proaching case is larger than the gain (loss) in the reced- 
ing case for an(a) attractive (repulsive) potential. The 
reason is that if we keep the setting the same but only 
change the direction of V, both v and n will change. In 
the attractive case (Figure^), if mi moves towards mi 
the angle between the asymptotes, 2ip, is smaller than 
it would be if mi moved away. Therefore, |n|= 2cos , i/> is 
larger in the approaching case. For very high speed en- 
counters, |vi| S> |va|, we can approximate v to be the 
same in the two cases, therefore Eq. [3] becomes: 



AKE mi = fivV.n 



-2fivVcost/j a 
2fjLvVcosip r 



(4) 

We next consider the case where mi moves in random 
directions. We do so by applying the above to the 



cases of gravitational and Coulomb interactions, where 
the equations of motion are particularly simple to 
integrate (Io| and analytical results can be compared 
to the numerical experiments of the previous section. 
We first parametrize: V2 = viCos9i + visin8j which 
is equal to V in the limit where mi 3> m\. The 
orbits of scatterings in V(r) = —a/r (here a is taken 
to be positive) are hyperbolas parametrized by Q: 
r = a(ecosh£ — 1), x = a(e — cosh£), y — a\/e 2 — lsinh£ 

and t = \J (esinh£ — £), where £ (hence time) 

varies from -co to oo and v = ff + We 
are interested in the change in the kinetic energy 
of the smaller mass in the lab frame [llj. Since, 
mi 3> mi =$-fi = mi; we have for AKE mi = 
lim t ^ 00 (KE mi ) - lim t ^_ oc (KE mi ), wh ich reads 
AKE mi = -IKy^pcosfl = -21/y^pcos7/.cos6>. 
Here 9 £ [0, 2ir) and 9 — corresponds to the approach- 
ing case. At first sight it seems like averaging over 9 
would yield zero but as mentioned above, the eccentricity 
e , the semi-axis a and consequently tp depend on the 
initial velocities as they are functions of the angular 
momentum and the energy in the center of mass frame. 

We can rewrite the formula by noting that costp = 1/e 
and a — mi ° 2 . For small angle scatterings, miDv^ 3> 

a, we have e w miD v?L, where Woo = wL. ^ is the 
relative velocity before collision and is a function of the 
initial angle the two velocities make with one another. 
Therefore the previous equation can be rewritten as 



AKE mi = -2V 



m\a. 
a 



cos tjj cos 9 



2aV 

DVoo 



cos 9, (5) 



with D being the impact parameter, which also de- 
pends on the initial angle between the velocity of the 
two masses. The strength of interaction for gravitational 
and Coulomb attractions are a = Gm\mi and a 
respectively. 

For the two cases discussed above we have for the ap- 
proaching case: 9 = and vq = v a . And for the receding 
case 9 = ir and vo = v r . It is easy to see that v a < v r ; 
therefore we yet again arrive at the same conclusion. 

In Figure 3] we plot Eq. O as a function of 9 and see 
that our theoretical treatment and numerical experiment 
Figure [2] agree. 

A formal statistical treatment depends on the particu- 
lar laws of interaction as well as the statistical distribu- 
tion of the large masses. For V = a/r it remains to show 
that under the assumption of microscopic chaos in the 
direction of motion of the large masses, the effect holds 
for any generic choice of initial conditions. Equation [5] 
can be written entirely in terms of initial conditions in 
the lab frame as (see Appendix IIII CI) 



AKE„ 



2aV 



2aV cos 9 
\p\ \sin6-q/p\ 



(0) 
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Theoretical Prediction: Energy vs. Theta 
x 1 (T 5 (True asymptotics are compared) 
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Figure 4: Theoretical prediction of a net loss of energy for 
attractive fields agrees with the numerical experiment. Here, 
was shifted by ^ to be comparable to the experiment shown 
in Figure [2] As before 9 = £ corresponds to the approaching 
case and 8 = to the receding case. The small discrepancy 
between the theoretical and experimental values of AE is due 
to the finite box size in the experiments. 



where, p = (x2 — x\) V2 and q = — (7/2 — Hi) Vi encode 
the dependence on the initial conditions. Now let z = q/p 
and (,) g denote the 9- average, we have 



(AKE mi ) g cx f(z) 



cos 9 



| sin 9 — z\ 



(7) 



where / encodes the dependence on the choice of initial 
conditions. Figure[5]establishes the existence of the effect 
irrespective of particular choices of initial conditions in 
the scattering process as long as the assumptions leading 
up to Eq. [5] are intact. 

One can use mechanics of similarity to calculate 
AKE mi when the law of interaction is Coulombic [8|. If 
we fix the corresponding paths we can related the change 
in the kinetic energy of the two: AKE c mi = SaAKEJ^, , 
where we have explicitly labeled the quantities corre- 
sponding to the gravitational and the Coulomb interac- 
tions by g and c respectively. Furthermore, the corre- 
sponding time of travel between two fixed points on the 

path are related by j- = y^-. Consequently, for repul- 
sive Coulomb interaction our final result reads, 



AKE mi = 



2a c V 



cos( 



(8) 



where a c = |^ as before. We see that the equation 
we get is the negative of Eq. [S] which, as expected, cor- 
responds to a net gain in energy when averaged over 9. 
A calculation from first principles yields the same result 
(see Eq[l3]in Appendix IB]). 

Statistical Context- A phe nomena worth considering 
is dynamical friction [lH. fl5|. Dynamical friction, how- 
ever, is like Brownian motion 12| as a big mass enters 
a medium of many smaller particles and slows down as 



Figure 5: Plot of / vs. z given by Eq. [7] 



a result. But it is distinct from Brownian motion as the 
interactions are long ranged. It is found that in dynam- 
ical friction "only stars with velocities less than the one 
under consideration contribute to the effect" [H| and @, 
p. 299]. 

The main requirement in our work is that the particle 
scatters from the time- dependent field of other massive 
particles, that are moving randomly and slowly in space, 
through a series of two-body scatterings. In addition, it is 
required that during the effective scattering the scatterer 
remains on one side of the scattering particle. This is 
equivalent to the assumption of the slow motion, in the 
dilute medium, of the large masses. 

It would be interesting to analyze the effect of micro- 
dynamics of the structure in the universe on the fre- 
quency shift of photons coming from distant sources. We 
expect a small loss of energy for photons that undergo 
weak lensing 1^, 2(3] in the interstellar media of galax- 
ies, clusters and super clusters [21]. Clearly the potential 
fluctuations, derived within this framework, in energy of 
the photons will be due to the dynamics of the lenses and 
disjoint from the essential time dependence of large scale 
space-time metric in the general theoryof relativity as in 
Sachs- Wolfe and Rees-Sciama effects (lMi- 

At first sight the effect we are proposing seems to vi- 
olate the equipartition of energy because a smaller and 
lighter particle "heats up" the medium of much larger 
particles that have higher average energies. But note 
that we are working with a non-equilibrium process in an 
open system In attractive potentials the small par- 
ticle starts from non-equilibrium initial conditions and 
through a series of scatterings it statistically loses en- 
ergy till a final scattering where the energy in the center 
of mass is negative. There on the small particle would 
have a bounded orbit about that final scatterer. This 
corresponds to the breakdown of small angle scattering 
assumption of this paper. In plasma physics this is known 
as shielding and in astrophysics it corresponds to captur- 
ing of a comet by a center of force. In repulsive potentials 
the energy of the small particle grows till the relativis- 
tic effects become significant and the transfer of energy 
between the small particle and the scatterers becomes of 
order unity with respect to the initial energy [22I The sec- 
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ond to last paragraph of section 13]. Hence the effect is 
not an equilibration process and is applicable to systems 
where the assumptions of small angle scattering, as well 
as, mi "C m 2 , V\ ^> v 2 but nevertheless m^j <C m 2 v\ 
hold. The former assumption is bound to break on time 
scales required for relaxation to equilibrium. 

It is our hope that this small effect helps with bet- 
ter understanding of rapid structure formation [23] , high 
energy cosmic rays, fermi acceleration |5j and redshift 
problem. 
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I. APPENDIX: NUMERICAL DETAILS 

The computer used is an IBM ThinkPad laptop. The 
experiments were done using an adaptive Bulirsch-Stoer 
algorithm [24] in C++. The simulation was done in three 
dimensions and I checked to make sure that the orbit 
stays in the plane and that the total energy is conserved; 
both to an accuracy of 0(1O -12 ). The data was then 
extracted and the plots were made in Matlab (Figure 
2). The theoretical plots were all generated in Matlab. 
As stated before, we investigated the case where — = 

' G 7712 

5 x ICR 8 , a = -0.13346, and k = 1. 

For numerical purposes it is convenient to rewrite the 
equation 



/ 777/1 Ct 

AKE = ~ 2Vw — — cos ip cos 9 (9) 

in terms of geometrical constants of the orbit, namely 
eccentricity e and a — Hence, e i— > e(9) ;a i— > a(9) 
and are given by the energy and the angular momentum, 
which are in turn function of the initial conditions. 
Therefore Eq. Oreads AKE mi = -^ v /^cos(9, 

where a = ^ and e = A /l + Here I and E de- 

note the magnitude of the angular momentum and the 
energy in the center of mass frame. 

A gravitational model for this choice of parameters 
is one in which (in MKS units) mi = 1.0, 7772 = 
2. x 10 9 , G = 6.67300 x 10~ n . Therefore we have 

= 7nim 2 = q.9999999995 and a = 0.13346. The 

' 7771 +7772 

big mass 7712 was initially at the origin and mi had 
the coordinates xq = —10000, i/q — —2000 with ini- 
tial velocity v x o = 0.008, v v q = 0.005. The initial ve- 
locity of m,2 was V2 = Vcos6>i + Vsin6>j; therefore, 
v o($) = (vxo — Vs'm9) i + (v y o — Vcos6>) j from which we 
can calculate the constants of the orbit. The angular 
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momentum is conserved and is given by, 1 = /zro x vo = 
/i [xo(v y o — ycos^) — yo(v X Q — V^sinfl)] k, and the energy 
is given by, E — i/i|v (#)| 2 — tj^t. Therefore the orbit 
constants a and e are specified and so is [SJ which is plot- 
ted in Figure 4 as a function of 9 and compared to the 
experimental case. 

The experiments were done by integrating Newton's 
equations of motions in the lab frame subject to the same 
initial conditions. Namely, I integrated 



ri 



T2 



Gm-2 

w 

Gmj 



0i - r 2 ) 



r 2 - ri 



(10) 
(11) 



using Bulirsch-Stoer integrator. I then saved the state 
of both particles (position and velocity) as well as the 
change in the energy of the whole system to make sure 
it was constant. By conservation of angular momentum 
motion must take place in a plane, which I took to be the 
xy plane and monitored the z component of the position 
of both bodies. Accuracies were 0(1O~ 12 ). Furthermore 
I calculated, in the C++ code, the kinetic energies of 
both particles for values of 9 = pj^, where p = 0, • • • , 19 . 
These were written in .txt files that in turn were uploaded 
in a Matlab code which plotted the asymptotic value of 
the kinetic energy of mi vs. 9. The result is Figure 2. 



AKE, 



2aV 



cos 9 



(13) 



III. APPENDIX: TRANSFORMATION 
BETWEEN THE LAB AND CENTER OF MASS 
FRAME 

A. Transformation from the Lab Frame to the 
Center of Mass Frame 

Given the initial conditions in the lab frame we wish 
to obtain the parameters of the orbit which fully specify 
the problem in the center of mass frame. 



In the Lab at to 


In the COM 


n = (xi,yi) 


E 


vi = (v x i,v y i) 


1 


*2 = (^2,2/2) 


D 


V 2 = (v x 2,V y2 ) 


Woo = |Vl - V 2 | i _, =too 



As before we have /x 



mi+m 2 



and at t = t : r = 



mi-t-m2 ' 

( r i — r 2) 5 v o = v i — v 2- The energy in the lab is given 
by E L = ^vi 2 + ^v 2 2 - U(r ) and is related to the 
energy in the center of mass by: 



II. APPENDIX: DERIVATIONS IN THE CASE 
OF COULOMB REPULSION FROM FIRST 
PRINCIPLES 



where, 



E = Et 



(mi + m 2 ) 



V? 



COM 



(14) 



The gain in the repulsive force field with k = 1 is done 
similarly. Let Q: 



r = a (ecosh£ + 1) , t — ^/ [ia? ja (esinh£ + £) 
x = a (cosh £ + e) , y = ciy/Je 



Consequently, v 



1) sinh£ 

a smh ^i+Ve 2 — 1 cosh 



and v 



fia ecosh£+l 
1 

M = (ffofMri) ' • As for the attractive case 

the motion of m 2 in the lab frame taken to be 
V = V+os0i+ V sm9']. Eq. [3] becomes: 



lim AKE mi = 2Vt — cos V 1 cos ( 

?->oo V a, 



cos-0 = 



(12) 



for small angle scatterings cos?/> ~ 0; therefore, cost)) 
a/miVgD <C 1. The final form becomes: 



V? 



COM 



m\V\ + m 2 v 2 
mi + m 2 
m\v\ + m 2 w| + 2mim 2 (vi.v 2 ) 



V 2 



(mi + m 2 y 
2/iviV2 cos 9 
(mi + m 2 ) 



(15) 
(16) 
(17) 



Vcom is used to emphasize that it is obtained from the 
initial conditions; otherwise it is equivalent to V. 

The last equality is obtained by assuming vi , with no 
loss of generality, to point along the x-axis; further Vq 
is independent of 9. Therefore, the energy in the COM 
reads, 



E = Eq — nv\V2 cos 9; where (18) 

(mi + ma) o 
E = E L - K 2 2 V 2 . (19) 

Angular Momentum: 1 = [ivq x vo denotes the angular 
momentum and relative velocity. At infinite separation, 
Woo, is determined by the energy through, E 



— T 



7 




Lastly, the impact parameter D is determined by the 
angular momentum and Vr*, via D = — — . With this, 
the problem is fully specified in the center of mass frame 
from the quantities given in the lab frame. 

B. From Center of Mass (COM) to The Lab Frame 

This is more standard and is done in most textbooks. 
See for example^ Section 14]. Though, what we said ap- 
plies to all central forces, recall that k = 1 corresponds 
to gravitational and Coulomb attraction. For k = 2 the 
integrals can be integrated in terms of elementary func- 
tions @, p. 40]. 



Note that the expression depends on the square of the 
angular momentum. More importantly note that the pre- 
factor to the cos 9 depends on 6. The angular momentum 
expressed in lab quantities reads, 

1 = mi [(xi - x 2 ) (-u 2 sin6>) - (y x - y 2 ) (vi - v 2 cos9)} 
~ mi [« 2 (x 2 - x\) sin - v x {y x - y 2 )} 

which when substituted into [2T] gives 



— 2Voy / 2miEo — m\v\v 2 cos 9 cos 9 



'mi — 




21 2 (Eq— mi v\ V2 cos 8) 
mia' 2 



If we further assume e 3> 1, the top equation Eq|2T1 
reduces to 



Formulation of Our Results Entirely in The Lab 
Frame 



Couple of remarks are in order before we proceed, 
first note that when m 2 ^> mi, fi = mi and therefore, 

V . 



VcOM = 



2m\v\V2 cos 6 
(mi+m 2 ) 

From section 8.3, we can obtain all the quantities in 
the lab frame, in particular we no loss of generality we 
take, vi = «ii, v 2 — v 2 cos^i + v 2 sin#j. This specifies 
v = V! - v 2 . In addition, r = (xi -x 2 )i+ (y x - y 2 )j- 
Hence, from Eqs. [S]and[S]we have, 



AKE, n 



2aV 



cos 9 



2aV 



cos 9 



(22) 



Dwoo " |p| |sin 6* — q/p\ 

where, p = (x 2 — Xi)v 2 and q = — (y 2 — y±) Vi encode 
the dependence on the initial conditions. Now let z = q/p 
and (.) e denote the 9- average. Therefore we have, 



(AKE mi ) e ^f(z) = 



sm — z\ 



(23) 



AKE„ 



-2am\Vo\ 



i (2Eq — miviv 2 cos 9) cos 2 9 
mia 2 + 2l 2 (Eq — miviv 2 cos 9) 



(21) 



here / determines the dependence of this effect on the 
choice of initial conditions. The plot of / vs. z is shown 
in Figure 5. 



